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The generation of the population oscillation of the multicomponent spinor Bose-Einstein conden- 
sate is demonstrated in this paper. We observe and examine the nonsynchronous decreasing pro- 
cesses of the magnetic fields generated between quadrupole coils and Ioffe coils during the switch-off 
of the quadrupole-Ioffe-configuration trap, which is considered to induce a nonadiabatic transition. 
Starting from the two-level Schrodinger equation, we have done some numerical fitting and derived 
an analytical expression identical to the results of E. Majorana and C. Zener, of which both the 
results well match the experimental data. 



PACS numbers: 32.80.pj, 03.75.Mn, 03.75.Kk 
I. INTRODUCTION 

Recently, the multicomponent spinor Bose-Einstein 
condensate (BEC) has become a "hot" topic, since it 
shows an appealing expectation in providing entangled 
spin systems applied in quantum optics and quantum 
computations Q, Q , and in the quantized votex applied 
in the study of superconductors and superfluidity 0, . 
Since more features [H IS @> @ an d more phenomena 
|ToL ITTl IrH Il3| are explored, the quantitatively manip- 
ulating and splitting the multicomponent spinor BEC is 
urgently desired. In several groups, the means of opti- 
cal manipulation has been applied to split and study the 
spinor BEC [14. Il5l ll6L Il7 | . However, since the BEC is 
more commonly generated in static magnetic traps, the 
splitting and manipulating by the means of magnetic field 
will be more convenient and efficient 10 . 

Nonadiabatic transitions of atoms within magnetic 
sub-levels are very old problems established in about 
1930s [H [H Hj. The model of Majorana transition 
in which the magnetic held evolves as B x (t) = 0, B y (t) = 
Const, B z (t) = kt and t = (—00,00) was created and 
solved by E. Majorana in 1932 [2]]]. Then, a more explicit 
and easily comprehensible explanation about the Majo- 
rana transition was described by I. I. Rabi j^. After 
that, though further theoretical discussions and analy- 
ses and other implements such as group theories have 
emerged to improve the study of Majorana transition 



@ HE El USUI, few quantitative experiments nor 
theoretical investigations have been undertaken, because 
there were not any precise ways to control the fleeting 
hot atoms, and some qualitative explanations about spin 
flips and atom loss in magnetic quadrupole traps are well 
known [29t I30L l3lj. However, with the development in 
the experiments and techniques of ultracold atoms, the 
almost motionless atoms can be provided to interact with 
the swiftly rotating magnetic field, and it makes possible 
the quantitative examination of Majorana transitions. In 
our previous work, we demonstrated that Majorana tran- 
sitions will emerge during the switch-off of the magnetic 
trap in our experiment system and induce the split mul- 
ticomponent spinor BEC [Tfj|. 

In this paper, we report the observation and the expla- 
nation of population oscillation of the multicomponent 
spinor BEC induced by nonadiabatic transitions. We be- 
lieve that the formation of the population oscillation is 
caused by the vertical oscillation of the condensate cloud 
and the nonsynchronous decreasing processes (NDP) of 
the magnetic fields, which are both measured experimen- 
tally and demonstrated in this paper. Further, we de- 
rive an analytical expression of this, starting from the 
Schrodinger equation and the experimental conditions, 
while the numerical calculation has also been undertaken. 
Both the analytical and numerical fitting results well fit 
the experiment data. 
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II. GENERATING THE POPULATION 
OSCILLATION OF MULTICOMPONENT 
SPINOR BEC 



Our experiment is set up on a standard equipment 
system for generating a cigar-like condensate in a di- 
lute gas of 87 Rb. A compact low-power quadrupole-Ioffe- 
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FIG. 1: (Color online) The Vertical Oscillation: The y- 
direction coordinate labels versus the vertical oscillating time 
in the trap. The vertical down is the direction of the gravity. 
Different oscillating frequency is due to different tightness of 
trapping confinement. The oscillation before the time 17.5 ms 
represents the switch-off of the compensate coils, while after 
it represents the restoring of the compensate coils. When the 
compensate coils are switched off, the average velocity of the 
oscillation is around I cm/s, and the period is calculated to 
be around 1 1 ms. 

configuration (QUIC) trap with trapping frequencies of 
u> r = 2tt x 220 Hz in radial directions and lo z — 2ttx 20 Hz 
in the axial direction, in which the quadrupole coils are 
assumed to be along x-direction and the Ioffe coils along 
the z-direction, is mounted on the lower chamber where 
the vacuum is up to 2 x 10 -11 mbar. After the evapo- 
rative cooling, the condensate cloud is held in the center 
of magnetic trapping potential generated by 23.6 A cur- 
rent in Ioffe coils and 24.1 A current in quadrupole coils. 
In addition, in order to reduce the offset field Bq in the 
center of the trap, we also apply a couple of compensate 
coils along the z-direction to regulate the tightness of the 
trap center. When the compensate coils are switched on, 
the Bq is reduced from 8.5 Gauss to 1.5 Gauss and the 
typical gradient of our trap is about 150 Gauss/cm. 

Firstly, the compensate coils are switched off, and af- 
ter a certain time interval, which we can call "coil delay 
time", the QUIC trap (the quadrupole coils and Ioffe 
coils) is switched off afterwards. During this coil delay 
time, as we know, the condensate cloud will be oscillat- 
ing up and down vertically due to the balance of gravity 
and new trapping potential whose confinement is much 
tighter than the original one due to the change of Bq. 
To demonstrate the vertical oscillation clearly, we first 
switch off the compensate coils and then restore it (see 
Fig. 1). We can see that, at different coil delay time dur- 
ing the oscillation, the condensate cloud may stay at the 
different altitude which can be expressed 

Y verUc ai(prn) = 23.2 x cos0.564i(ms) - 28.2, (1) 

where the positive direction of y-direction is vertically up 
(opposite to that of gravity) . From both the expression 
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FIG. 2: (Color online) The Measurement of Nonsynchronous 
Decreasing Processes (NDP) : Induced voltage of the detective 
coil when the quadrupole and Ioffe coils are shut down versus 
time during the decreasing process of the currents. We applied 
a detective coil to examine the evolution of the field generated 
by the two types of coils respectively. While one is measured, 
the other is substituted by an inductance vessel. The blue 
line represents the detecting results of quadrople coils, while 
the red one represents that of Ioffe coils. The time constant 
of exponential evolution is respectively about 91 ms and 39 
ms. 



Eq. 0} and the Fig. 1, we find that the average veloc- 
ity of the oscillation is around 1 cm/s and the period is 
calculated to be around 11 ms. 

The fact is, though the quadrupole trap and Ioffe trap 
are switched off at the same the time, their decreasing 
processes are totally nonsynchronous according to our 
experimental measurement (see Fig. 2). We apply a de- 
tective coil to detect the field evolution of the two types 
of coils respectively. While one is measured by the detec- 
tive coil, the other is substituted by an inductance vessel. 
This means of detection can not give the time constant 
of evolution of the real field but the changing trend of 
the field. Due to the derivative feature of the exponen- 
tial function, the time constant we detected can be seen 
to be identical to that of the real field. 

As our previous work shows ^3], when the NDP of 
magnetic fields emerges, the condensate atoms will ex- 
perience the zero field B z = and the reversion of the 
direction. In this case, there comes a nonadiabatic transi- 
tion among magnetic sublevels, and the different compo- 
nents of spinor BEC will be separated in space by Stern- 
Gerlach Effect due to the interaction between atom spins 
and the gradient of magnetic field. The population os- 
cillation is shown partly from 18ms of coil delay time to 
28ms of coil delay time (see Fig. 3), and surely much 
longer oscillation has been observed in our experiment. 
Examining all the pictures during the whole process, we 
find out that the period of the population oscillation is 
about 11 ms. The total number of the atoms is about 



3 



2 x 10 5 . The most left cloud is considered to represent 
the vnp = +2 component and the most right one belongs 
to mp = —2 component, for we assume that the direc- 
tion of the field after the reversion is the quantum axis. 
It is necessary to point out that the periods of the ver- 
tical oscillations and population oscillation are more or 
less the same. We believe this means that the population 
oscillation is induced by the vertical oscillation and the 
different distributions of atoms are related to the differ- 
ent altitudes of the condensate cloud. 



III. THE BASIC MODEL FOR TWO-LEVEL 
NONADIABATIC TRANSITION 

According to I. I. Rabi's description, the Majorana 
transitions only take place in the case that the rotat- 
ing frequency of the magnetic field fiiot. = dB(t)/2irdt 
is big enough to be comparable to the Larmor frequency 
of the field /i ar . = g/ioB(t) /2ttH \22\. This is easy to un- 
derstand from the perspective of magnetic resonant tran- 
sitions (MKT). Because the emergence of the zero field 
and the reversion of the direction often bring the small 
magnitude of the field and huge rotating frequency, the 
nonadiabatic transition will happen more easily. In our 
experiment, the evolution of the magnetic field just co- 
incides with this conclusion (see Fig. 4). When the field 
reverses its direction, the rotating frequency is remark- 
ably huge. At this time, the magnetic moment of the 
atom cannot follow the rotating field and there will be 
a transition among the magnetic sublevels that can be 
ascribed to the Majorana transition. 

Majorana Formula has been derived from both quan- 
tum mechanics and group theories to explain the multi- 
level cases [^1 128|. For a system with a total angular 
moment J, 



Pm,m> = (J + m)! (J + to')!(J - m)!( J - m')\ \ cos - 
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FIG. 3: (Color online) The Population Oscillation: Images 
of multicomponent spinor BEC due to different "coil delay 
times" . Each component of spinor BEC is separated by Stern- 
Galach effect due to the interaction between atom spin and 
the gradient of the magnetic field. The bar charts of popula- 
tion distribution of each image are also shown. From a) to f), 
the pictures are one period of the oscillation and the period 
value is examined to be around 11ms. The color bar shows 
the optical density of each pixel on the absorption images. 
We consider the most left cloud representing the mp = +2 
component and the most right one representing the rriF — —2 
component, for we assume the direction of the field that is 
after the reversion to be the quantum axis. 



where the value of the parameter 9 is given by the two- 
level transition 



sin" - = Pi 



(3) 



This means that once we solve the two-level case, the 
results can be generalized for any system with any value 
of J. So, in this part, we will focus on the physical model 
of the two-level nonadiabatic transition. 

The time, during which fft, t. is big enough to be com- 
parable to the jhar.i is about 1 /is (see Fig. 4a), so the 
movement of the cloud is about 0.01 /xm according to the 
average velocity about 1 cm/s. Compared with the ver- 
tical oscillation magnitude (~ 50 /Ltm), atom cloud can 
be seen to be motionless. So we consider a system of 



motionless atoms whose spin moment is s = ~ with a 

time evolving magnetic field B(t). In this simple case, 
we begin with the Schrodinger equation 



H \ C2 



(4) 



where we know F = |<? in the two-level case. According 
to the configuration of QUIC trap, the center of the mag- 
netic trap is right on the z-axis. In addition to the drag- 
ging down of the gravity, however, the center of the to- 
tal trapping potential is slightly down along y-direction. 
So, the three components of the evolving magnetic field 
should take the form B{t) = (0, B y (t), B z (t)). Putting 
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FIG. 4: (Color online) a) The Comparison of The Two Fre- 
quencies: The Larmor frequency fhar. is proportional to the 
magnitude of the field, so when the field reverse its direction, 
it is small. Meanwhile, the rotating frequency /.Rot. is re- 
markably huge due to the reversion. The time, during which 
fnot. is big enough to be comparable to the fhar., is about 1 
fj,s. b) The Real Evolution of The Magnetic Field B z : When 
the QUIC trap is switched off, the zero field B z = emerges 
around 5 /is. Though the decreasing processes of the field 
are exponential due to the discharge of the coils, the part at 
which the transition happens is approximately linear. From 
these two figures, it is confirmed that the emergence of Bo = 
and the reversion of the direction will induce the nonadiabatic 
Majorana transition. 



the matrix form into the basic equation Eq. J3J and tak- 
ing the substitution a = g/is/^h, one yields 



c'i = —iaB z ci — aB y C2, 
c 2 = aB y ci + iaB z c 2 , 



and the initial conditions are: 

Cl (0) = 
c 2 (0) = 



(5) 



(0) 



Substituting the real evolution of the magnetic field into 
Eq. JSJ, one can get the transition probability. 

Though, as we know, the evolution of the magnetic 
field of the coils is exponentially down due to the dis- 



f B y (t) = B yi exp {—nt) + B yq exp (— r q t) . 
\ B z (t) = B zi exp (-nt) - B zq exp (-r q t) , 



(7) 



where B y i and B Z i are respectively the y-direction and 
z-direction component of magnetic field generated from 
the Ioffe coils, B yq and B zq are from the quadrupole 
coils, and Tj and r q are the reciprocals of the exponen- 
tial time constants Tj a ff e and TQ ua d in Fig. 2. Since 
only the transition properties will be taken into account, 
it is appropriate to take the first order approximation 
exp (— T Ql it) ~ 1 — T q ^t and describe them in the linear 
forms (see Fig. 4b) 



By{t) 



Ay - Cyt, 

A z — C z t, 



where 



A z = (B zl 

Cy 

c z 



(B zi Ti 



Byq ) , 

Bzq) , 
+ B Vq T q) 
- B z q T q ) 



(8) 



(9) 



From above we know Tj « l/39ms is twice of r q 
1/9 lms (see Fig. 2) and hence tj 



2r q . So the value of 



C y and C z 
and C z « i 
satisfy C z 



can be estimated to be C v m r a (2B V 



Byq) 



■ q (2B zi ~ B zq ), and the typical values of theirs 
^> Cy. So we can simplify the expressions 



of magnetic fields by fixing the B y at the value A y $ = 
A y — Cy ■ to at the time t = to when the field B z reverses 
its direction at B z (t Q ) — 0, 



By(t) 
Bz(t) 



A z - 



C z t. 



(10) 



With Eq. (|l(Jf> and all other substitutions taken into the 
Eq. JSJl, we can get the second-order differential equa- 
tions which can be transformed into Webber equations 
[2fj| . As referred above, the transition only occurs when 
B z reverses its direction, so it is allowed to reset the ini- 
tial condition so as to utilize the asymptotic solutions of 
Webber equations 



ci(-oo) = 1, 
c 2 (-oo) = 0. 



(11) 



At the same time, the transition probability should take 
the form P = 1 — |c2(+oo)[^ since the magnetic field has 
reversed its direction H3,l^| . 

Therefore, we can derive the analytical expression of 
the Webber equation with its infinite asymptotic solu- 
tions. If we take this substitution 



aAl 



2C Z 



(12) 



the final result of the analytical expression is (the calcu- 
lation details can be found in |20| and [32| ) : 



Pi 



exp(— 27ra). 



(13) 
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FIG. 5: (Color online) The Experimental Data and Analyzing 
Results: The population oscillation of \F = 2, ulf = +2} state 
versus coil delay time. The red circles represent the experi- 
mental data, while the blue triangles represent the numerical 
fitting. The dashed line is drawn by the analytic solution from 
above theoretical deducing. 



The above analytic expression is identical to the results 
derived previously 20,21]. From this analytical solution, 
the population oscillation should rely on the magnitude 
of A y o which is the only unfixed parameter in a. In our 
experiment, A y o is approximately proportional to alti- 
tude oscillation Y ver ucai , because the value of B V i + B yq 



is around linear in radial directions. Hence, we approxi- 
mately have a w kY^ ertical and A: is a constant. Finally, 
we know that the transition probability will also oscillate 
with the altitude oscillation Y vert i ca i 



IV. CONCLUSION 

To conclude, we observed the split spinor BEC and 
the population oscillation of the five components of 87 Rb 
generated in the magnetic trap. For the first time, ac- 
cording to our knowledge, the population oscillation in 
which the population distribution of the five components 
can be regulated by setting the experiment parameters 
is reported and explained for its experimental formation. 
We believe that our experiment could most probably offer 
a novel and profound means to produce and manipulate 
the spinor BEC and the distributions of its components. 
In addition, we analyze the phenomenon of BEC split- 
ting, which coincides with the explanation of NDP pro- 
posed by ANU group j3^|, and measure the value of the 
DNP experimentally. 

Starting with Schrodinger equation, together with the 
experimental conditions, we derive the analytical solu- 
tion of nonadiabatic Majorana transition of condensate 
atoms within the magnetic sublevels. Though Majorana 
transition has been broadly applied in cold atom physics 
for qualitative explanations, it is for the first time that 
we quantitatively apply it to well fit our experimental 
data of the nonadiabatic transition. 

For further experiments, we have designed a more com- 
plex system to control the magnetic fields of the two 
types of coils separately, and with the aids of the the- 
oretical analysis mentioned above, the arbitrary manipu- 
lations of the atom populations among the multicompo- 
nent spinor BEC are anticipated. Besides, we also plan 
to load the separated multicomponent condensates into 
the optical trap and study the interaction between the 
different components. To summarize, our experimental 
method provides a convincing and potential way of study- 
ing the ultracold atom physics, including BEC and atom 
optics. 



Pi 



exp (- 



-2irkYv ertical ) 



(14) 
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